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We propose a model for the pair of electromagnon excitations observed in the class of multiferroic 
materials 7?Mn03 (R is a rare-earth ion). The model is based on a harmonic cycloid ground state 
interacting with a zone-edge magnon and its twin excitation separated in momentum space by 
two times the cycloid wave vector. The pair of electromagnons is activated by cross coupling 
between magnetostriction and spin-orbit interactions. Remarkably, the spectral weight of the twin 
electromagnon is directly related to the presence of a magnetically induced oscillatory polarization in 
the ground state. This leads to the surprising prediction that TbMnOa has an oscillatory polarization 
with amplitude 50 times larger than its uniform polarization. 

PACS numbers: 75.80.+q, 78.20.Ls, 71.70.Ej, 75.30.Et 



I. INTRODUCTION 

The coexistence of magnetic and ferroelectric phases in 
multiferroic materials gives rise to important effects re- 
lated to cross correlation between order parameters and 
external fields*^ 2 - A notable consequence is that the el- 
ementary excitations are not purely magnetic or ferro- 
electric in character. The spin waves are admixed with 
vibrational modes of the electric polarization, giving rise 
to electric dipole active magnons, the so-called electro- 
magnons. Electromagnons were postulated to exist in 
1969 by Bary'akhtar and Chupis 3 but their existence 
was confirmed only recently by sensitive optical experi- 
ments in the far infrared^^ 7 -!^ combined with neutron 
scattering ^ 

Apart from its fundamental importance, the electro- 
magnon spectra gives invaluable information on how 
magnetism couples to ferroelectricity. The class of per- 
ovskite manganites i?MnC>3 (RMO) is playing a major 
role in this respect. Here R is a rare-earth ion such as 
Gd, Tb, Dy, or mixtures between them (e.g., Gd^Tbi-^, 
etc.). At low temperatures (T ~ 20 K) the Mn spins 
are ordered in a spiral state with period incommensu- 
rate with the latticed The spin spiral is of the cycloid 
type. This breaks space inversion and gives rise to a 
ferroelectric moment along one of the directions in the 
cycloid planed At the microscopic level, this ferroelec- 
tric moment originates from Dzyaloshinskii-Moriya (DM) 
coupling^i^^ 

Nevertheless, contrary to predictions based on DM 
coupling^ electromagnons in RMO are observed only 
when the electric field of light is along the crystallo- 
graphic direction ajZiSiSiil Moreover, for all ions R, two 
electromagnons are always observed^ Recently, the ori- 
gin of one of these electromagnons was explained by 
Aguilar et al.r& who pointed out that the high fre- 
quency electromagnon was a zone-edge magnon activated 
by pure magnetostriction. However, the origin of the low- 
frequency electromagnon remains unknown. 

A recent theory-^ shows that the anharmonicity of a 
distorted cycloid ground state gives rise to multiple elec- 



tromagnons that were observed using Raman scattering 
in BiFeOs. 6 However, the anharmonicity mechanism can- 
not explain the pair of electromagnons in RMO because 
the degree of anharmonicity detected by neutron diffrac- 
tion is too low in these materials^ 

Here we propose a model that is able to explain the 
double electromagnon feature observed in RMO. We 
show that a simple harmonic cycloid ground state with 
wave vector Q has a special pair of twin electromagnon 
excitations, located at the zone edge and 2Q away from it. 
Remarkably, the activation of the latter electromagnon is 
directly related to the presence of an oscillatory polariza- 
tion with wave vector 2Q in the ground state. 

II. MODEL 

The usual phenomenological Landau theory is based on 
a free-energy expansion into powers of the spatial deriva- 
tives of the order parameters^ Therefore, the Landau 
approac h 3 i 12 i 19 i 21 cannot describe large wave vector ex- 
citations such as zone-edge magnons. Here we adopt a 
model that is more microscopic than the Landau theory. 
Our model Hamiltonian consists of three distinct contri- 
butions, H — Hs + -ffph + H mo . The first contribution 
describes spin frustration in RMO^ ' 10 ' 11 ' 22 

H S =Y,Jn,rnS n -S m +DY / (Sn-a) 2 . (1) 

n,m n 

Here J n ,m is the exchange coupling between Mn spins at 
positions R n and R m , and D is a single- ion anisotropy 
along the crystallographic direction a (we assume D > 0, 
i.e., easy-plane anisotropy). RMO has an orthorhombic 
lattice, with an orthogonal system of crystal axis a, b, c 
and unequal bond lengths a =/= b ^ c . There are four 
Mn spins per unit cell; spins 1 and 2 lie in the ab planes, 
with spins 3 and 4 a distance c/2 above them. The cou- 
pling between nearest neighbor spins in the ab planes is 
ferromagnetic, and will be denoted by Jo; the coupling 
between ab layers along the c direction is antiferromag- 
netic and is denoted by J c . Spin frustration arises due 
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FIG. 1: (Color online) Coupled spin waves and optical 
phonons in the cycloid ground state. In-plane spin fluctu- 
ations are denoted by <f>, while out-of-plane spin fluctuations 
are denoted by ip. 

to coupling between next-nearest neighbors along the b 
direction; this is antiferromagnetic and denoted by J 2 b- 
When the stability condition J 2 b > — Jo I 2 is satisfied, 
the ground state is a be cycloi d 13 ' 14 i 15 as observed in 
experiments,—^ 

5o(R) = ±S"[cos(Q • R)6 + sin(Q • R)c]. (2) 

Here Q = Qb is the cycloid wave vector, with 
cos(Qb/2) = — J /(2J 2 fc). The upper (+) sign applies 
to ab layer spins 1 and 2, whereas the lower (-) sign 
applies to spins 3 and 4 in the neighboring ab layers im- 
mediately above and below them. Figure [1] depicts the 
cycloid ground state. The second ingredient of our model 
is the lattice fluctuation Hamiltonian, 

H P h = \m* (*» + w o*n) -e*Y, x n-E, (3) 

n n 

where relative displacement between cations and 

anions within each unit cell. Here e* is the Born ef- 
fective charge, m* is an effective mass, and loq is the 
(bare) phonon frequency. The modes x n are directly 
related to the polarization order parameter in Fourier 
space, P q (u) = -^fdt^e^-^-^Xn, where v 
is the volume of the unit cell. The electric field of 
light, Ee~'" Ajt , couples linearly to the q = polarization, 
Pq{uj) = Xoi^E, with electric susceptibility Xo( [j -') = 
{e* 2 )/[{ m *v,){u 2 -u 2 )]. 

The third contribution H mc models the coupling be- 
tween the spin degrees of freedom and the lattice. Fol- 
lowing the experimental observation that electromagnons 
are only excited when the electric field is along a£&iLll 
we consider the interaction, 

Hmc — e x n\9c{S l n — S ljn+b ){S 2n + S 2n+a ) 

n 

+9b(S\ n — S\ n+b )(S\ n + <§2,n+a) 

+(1-3,2-4)], (4) 

which is invariant under the Pbnm space-group opera- 
tions of RMO. Here g c ^ g b are coupling constants with 
dimension of electric fields; we neglect the components 



proportional to g a because these play no role in the dis- 
cussion below. Remarkably, Eq. @ is not rotationally 
invariant in spin operators. Therefore it differs in an 
important way from the pure magnetostrictive coupling 
considered previously^ The absence of rotational invari- 
ance in magnetoelectric coupling is a direct consequence 
of cross coupling between magnetostriction and spin-orbit 
interactions, and should be common property of all multi- 
ferroics with orthorhombic lattice (this follows from sym- 
metry since a ^ b ^ c precludes rotations that take one 
axis into the other). Later we will explain why the ex- 
perimental observations favor this coupling over other 
symmetry-allowed explanations. 

It turns out that Eq. ^ gives rise to a remarkable 
static effect. Minimizing i? p h with respect to x n , and 
plugging in the cycloidal spin order we get 

LEU = p IOP S in [(2n + l)Qb]a. (5) 

v o 

Hence the local polarization per unit cell is oscillatory 
with wave vector 2Q; the amplitude for the incommen- 
surate oscillatory polarization (IOP) is given by 

Piop = 4 (g b -g c ) X oS 2 sin (Q^J. (6) 

Remarkably, the presence of the IOP is directly related 
to the lack of rotational invariance in Eq. (|4|), i.e., the 
fact that gb ^ge- 
lt is well known that TbMnC>3 has a magnetically- 
induced lattice modulation with wave vector 2Q in the 
cycloidal phase. - However, we are not aware of any claims 
that such a modulation leads to oscillatory electric polar- 
ization. The main result of this paper is to show that this 
oscillatory polarization is directly related to the spectral 
weight of the twin electromagnon; hence it can be de- 
tected by optical experiments. 



III. TWIN ELECTROMAGNON EXCITATIONS 

We now consider the elementary excitations of the cou- 
pled spin-phonon system. The spin excitations 5S = 
S — So are parametrized as follows, 

SS i<n = tp iin a±^i tn [cos (Q ■ Ri in )c— sin (Q ■ Ri, n )b], (7) 

where the sign convention is the same as in Eq. Here 
the operators ip describe spin fluctuations out of the cy- 
cloidal plane, and <j> denotes in-plane (tangential) fluctu- 
ations (Fig. [T]) . We carry out a mean-field expansion of 
the Hamiltonian H about its equilibrium value by keep- 
ing only terms quadratic in the fluctuation operators, 
e.g., SP 2 , cj) 2 , SPtp, etc. Using the canonical commu- 
tation relations [</>j,n>i>k,m] = i^jk^nmS, we derive the 
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FIG. 2: (Color online) Typical dispersion curves along the 
cycloid direction b for the cyclon (blue solid) and extra-cyclon 
(red dashed) modes in a be cycloid. Our theory gives rise to 
twin electromagnons depicted by a circle and a square. When 
the magnetoelectric interaction is rotationally invariant [gi, = 
g c in Eq. @] , only the one denoted by a square is excited; in 
the absence of rotational invariance, both electromagnons are 
excited. Here Q is the magnitude of the cyclon wave vector 
and b the lattice constant along b axis. 



coupled equations of motion for spin and polarization, 



- fi C,q) {<t>U + fcq + 03<? + <K) = T g-k„ , (8a) 



- nt 



C,q+k Q J {<hq ~ foq + foq ~ 4>iq) = T q , (8b) 
(u 2 ~ n| C)q ) (<t>lq + <t>2q ~ 03? ~ <M = 0, (8c) 
) (01, " 02 9 ~ foq + <j>4q) = 0. (8d) 



•A/ 
^EC : q+k 0/ 



Here 



are Fourier transforms of the fields 
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iden- 



tical equations hold for the fields ipi n ^ There are four 
spin-wave modes, each transforming according to four 
different Id representations of the Pbnm space-group. 
The second mode [Eq. ([8b])] is related to the first mode 
by q — » q + ko, with ko = (2tt /b)b the Brillouin zone 
edge (the lattice has periodicity 6/2 along b). This corre- 
sponds to the fact that "anti-phase" fluctuations of the 
neighboring spins with the wave vector q are equivalent 
to "in-phase" fluctuations at q + 2n/b. Their dispersion 
satisfies fic,q->o — > 0, reflecting the phase sliding symme- 
try of the cycloid. Adding a constant phase (p to Eq. @ 
yields a different ground state with the same energy. At 
q = the first mode is a pure phase fluctuation (<pi n is 
the same for all i and ipi n = 0). For this reason this 
mode will be referred as a cyclon, with dispersion fic,q- 
The second mode is naturally called a zone- edge cyclon. 
Similarly, the fourth mode can be obtained from the third 
by the same translation in q space, q — > q + ko. However, 
the third and fourth modes are gapped (gap proportional 
to J c ). These modes will be denoted extra cyclons. The 
(bare) magnon dispersions are shown in Fig. [21 For q 



a) q = k 



b) q = k -2Q 





FIG. 3: Diagrams for twin electromagnon excitation by light, 
involving momentum exchange with the cycloidal ground 
state, (a) Zone-edge electromagnon at q — ko and (b) its 
twin at q = ko — 2Q. The fact that momentum is not con- 
served by fco = 2ir/b reflects the underlying lattice symmetry, 
which produces spin- wave modes connected by q — ► q + ko- 
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Typical values for the exchange couplings are J ~ —0.4 
meV, J 2b ~ 0.3 - 0.4 meV, and J c - 0.5 - 1 meV4 In- 
terestingly, only the zone-edge cyclon is coupled to the 
polarization. This occurs through the dynamical magne- 
toelectric field, 



-—7rb sin — — 



Qb 



h2 . D sm ^_^ 2 Jo^cos(J 

J 2b [cos (qb) - cos (Qb)] - (D 
{2(g b + g c )P q a (u J ) 

(9b -gc) [P q a - 2Q H + Pq+2Q^)}} 



Qb 
2 



2J C )} 



(10) 



Similarly, the equation of motion for the polarization 
only couples Pq to the zone-edge cyclon at q = and 
q = ±2Q. Hence we have a set of twin electromagnons 
at wave vectors q = ko and q = ko — 2Q. For g b = g c , only 
one electromagnon (the zone-edge cyclon) is activated^ 
This excitation corresponds to a scattering process with 
in and out-going momenta equal to the cycloid wave vec- 
tor Q [Fig. (3fa)]. However, when g b ^ g c , the twin elec- 
tromagnon (a cyclon at q = ko—2Q) appears, correspond- 
ing to two outgoing momenta adding up to 2Q [Fig. EJb)] . 



IV. DIELECTRIC FUNCTION AND 
COMPARISON TO EXPERIMENTS 

In order to relate to optical experiments, we calculate 
the dielectric function e(ui) — + AttPq (uj)/E a and 
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( s A 



Exp 

Th 



Gd .5Tbo. 5 Gdo.iTbo.g Tb Tbp.4tDyo.59 Dy 

0.30 0.39 0.40 0.36 0.37 

0.55 0.53 0.52 0.52 0.54 

0.08 0.11 0.13 0.41 0.51 



L(Si,-9c) 2 S C:ko J Th 


0.14 


0.16 


0.16 


0.16 


0.18 













TABLE I: This table compares our theoretical calculations to experiments in i?Mn03. A and B are, respectively, the lower- 
and higher-energy electromagnon modes detected in 



study its resonances. When the resonance frequencies 
are not too close to each other e is a sum of Lorentzians, 
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The electromagnon frequencies are seen to be down- 
shifted from the bare magnon frequency fic,q to 

— A 2 , while the phonon frequency ujq gets up- 



shifted to y <x>o + ^ph- The spectral weights S q share a 
simple relationship with the frequency shifts: For clcc- 



tromagnons, S q — [e(0) — EoJA^, while for the phonon 
<S P h = [e(0) — £00] i^Q — Ap h ) . The frequency shifts are re- 
lated by A* h = A* o + A| o _ _ 2 q, ensuring the satisfaction 
of the oscillator strength sum rule, Sk + <Sfc -2Q + <?ph = 
WqXo- The electromagnon spectral weights are given by 



J 2b 



(12a) 



o _ 27rSxo^o"o(gt>-gc) 2 tan 2 ^ / ^c,k Q -2Q 

fc °" 2Q ~ K-O 2 , )t0 _ 2Q ) 2 [cos^+-n 4 ^J I J2 " 



(12b) 



An important experimental result is that the spec- 
tral weights are nearly the same for the ab and be 
cycloids! 7 ' 9 ' 17 It is easy to see that this result follows 
from our model when \gb\ ^> |g a |,|<?c|- Note that in this 
limit the magnetoelectric coupling Eq. ([4]) is invariant 
under the flip of the cycloid plane, in spite of the fact 
that it lacks rotational symmetry. Table I compares our 
theoretical calculations to the measured values £Ia/^b 
and Sa/Sb, where A and B label the lower- and higher- 
energy electromagnons observed in RM04 We used the 
same model parameters as the ones in Fig. 4(b) of Ref. @ 
(filled symbols). The variations in the ratio Sa/Sb in- 
dicate differences in the couplings g c ,gb for different ions 
R. For Tb our theory matches the observed values when 
g c /g b = 0.05; for Dy we get g c /g b = -0.25. 



V. DISCUSSION AND CONCLUSIONS 

We now consider other possibilities for the activation of 
the low-frequency electromagnon. Interactions contain- 
ing crossed terms such as x^SfS^ are also allowed by the 
Pbnm symmetry. For a be cycloid, this interaction gives 
rise to an electromagnon at q = ko — 2Q. However, when 
the cycloid is flipped to the ab plane, this interaction 
leads instead to an electromagnon at q = ko — Q, with no 
electromagnon at q — k$ — 2Q. This result implies a large 



shift in electromagnon frequency (more than 50% as seen 
in Fig. [2]), that is in contradiction to experiments! 7 -^^ 
Therefore the interaction x^S^S^j can not explain the ori- 
gin of the low-frequency electromagnon. Similarly, in- 
teractions of the form x^SfSj only give rise to electro- 
magnons at ko — Q; this leads to ratios of electromagnon 
frequencies Qa/^b that are twice as large as obtained 
experimentally (Table 1)4 Hence x^SfSj is also ruled 
out. 

Aguilar et alX suggested that an elliptical spiral struc- 
ture such as the one in BiFeC>3 (Refs. y andKL9l) may ex- 
plain the low-frequency electromagnon. This possibility 
is also ruled out because the degree of ellipticity observed 
in neutron-scattering experiments^ is too low to explain 
the large spectral weight of the low frequency electro- 
magnon. Furthermore, Ref. suggested that purely mag- 
netic interactions such as SfSj would mix extra-cyclon 
magnons at q — Q to the zone-edge electromagnon, pro- 
viding an alternative explanation for the low frequency 
resonance. However, the extra-cyclon magnon energy at 
q = Q is approximately constant for different i?'s, in 
contradiction to the trend observed in experiments [Qa 
tends to decrease with increasing ionic radius; Qa for 
Gd . 5 Tbo.5 is 40% higher than for Dy (Ref. i)]. Inter- 
estingly, our theory gives a natural explanation for this 
trend since the value of the cyclon energy at q = ko — 2Q 
also decreases appreciably as the ionic radius increases. 
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We now discuss the implications of our model for the 
understanding of multiferroic order. One remarkable con- 
sequence of the presence of the q — k — 2Q electro- 
magnon is that its spectral weight can be directly re- 
lated to the presence of the IOP in the ground state. Us- 
ing Eqs. © and (fT2b|) we derive an important relation 
between the amplitude of the IOP and the twin electro- 
magnon spectral weight, 



p2 
MOP 



85 3 J- 



26 2 
COS 



TTVq 



Qb 
2 



Qb 



Qb 



S k 



-2Q 



kn 



-2Q 

(13) 

Hence optical experiments are a direct probe of the IOP 
amplitude. For TbMn03, the measured twin electro- 
magnon frequency is f2c,k -2Q — 25 cm -1 , with spectral 
weight S ko -2Q = 1-7 x 10 3 cm -2 . 8 Using Eq. flT3]) we ob- 
tain Piop = 4 /iC/cm 2 , a value that is 50 times larger 
than the uniform polarization Pq — 8 x 1CU 2 /iC/cm 2 
present in TbMn03. 

In conclusion, we introduced the concept of the twin 



electromagnon in order to explain optical experiments 
in the i?Mn03 family of multiferroics. Our symmetry 
analysis shows that an incommensurate oscillatory po- 
larization coexists with the well-known cycloid phase in 
these materials. Remarkably, there exists a direct re- 
lation between the twin electromagnon spectral weight 
and the amplitude of this oscillatory polarization. Hence 
we showed that TbMn03 has an oscillatory polarization 
with amplitude 50 times larger than its uniform polar- 
ization. This surprising conclusion underlines the impor- 
tance of electromagnons in the characterization of multi- 
ferroic order. 
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